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In this talk we will discuss the physics reach of the atmospheric neutrino data collected by 
a future megaton-class neutrino detector. After a general discussion of the potentialities 
of atmospheric neutrinos on general basis, we will consider concrete experimental setups and 
show that synergic effects exist between atmospheric and long-baseline neutrino data. Finally, 
we will show that present Super-Kamiokande data already have the capability to allow for a 
direct and unbiased measurement of the energy spectrum of the atmospheric neutrino fluxes. 



1 Introduction 

Despite their pioneering contribution to the discovery of neutrino oscillations, it is in general 
assumed that atmospheric neutrino data will no longer play an active role in neutrino physics in 
the coming years. This is mainly due to the large theoretical uncertainties arising from the poor 
knowledge of the atmospheric neutrino fluxes, which strongly contrast with the requirement of 
"precision" needed to further enhance our knowledge of the neutrino mass matrix. In this talk, 
we will show that despite these large uncertainties the atmospheric neutrino data collected by a 
megaton-class detector will still provide very useful information. 

Let us begin by reviewing what we have learned so far about neutrino oscillations. The solar 
and atmospheric mass-squared differences are clearly determined, and we know that the solar 
angle is large but non-maximal while that the atmospheric angle is practically maximal. The 
present best-fit point and la {3a) ranges areP 

^12 =33.7 ±1.3 {tti) , Ami =7.9^°:^ {tli,) x IQ-' e\\ 

O23 = 43.3 {+1-1) , |Am|i| = 2.6 ± 0.2 (0.6) x lO^^ eV^ , (1) 

^13 = 0^^:2 {tl'o') , 5cpG [0,360], 

leading to the following values for the elements of the leptonic mixing matrix, U, at 90% CL: 

0.53 ^ 0.58 

|[/|go^^ = I 0.32 ^ 0.49 0.52 ^ 0.69 0.60 ^ 0.76 | , (2) 

0.47 ^ 0.64 





and at the 3a level: 



0.50 0.61 0.00 ^ 0.20> 
|J7|3<7 = I 0.25 ^ 0.53 0.47 ^ 0.73 0.56 ^ 0.79 | . (3) 

0.42 ^ 0.69 0.61 ^ 0.83> 
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Figure 1: Allowed regions at 95% CL from the analysis of different combinations of neutrino experiments. In 
particular, the gray regions correspond the combination of solar, reactor and accelerator data. 

To set the basis for the following discussion, it is useful to study how much the information 
coming from atmospheric neutrino data contributes to this picture. The "solar" parameters 
9i2 and Arn^i are completely determined by the solar and KamLAND data alone, and will 
therefore not be discussed here. As for the remaining parameters, in Fig. [T] we show the allowed 
regions implied by different combinations of neutrino experiments. It is clear from this figure 
that the bound on ^13 comes mainly from Chooz, with a small contribution from solar and 
KamLAND data, while the atmospheric mass-squared difference Amgj^ is mainly determined by 
the accelerator experiments K2K and Minos. The only parameter whose determination is still 
dominated by atmospheric data is 623, being mainly a matter of total statistics. So it seems 
that indeed atmospheric experiments already have not much to say. 

However, at a better look we note that present reactor and accelerator data (gray regions) 
exhibit a very high degree of symmetry. In particular, they have practically no dependence on 
6cp, and they are totally insensitive to the neutrino mass hierarchy (sign of Am'^i) and to the 
octant (sign of 623 — 45°). Conversely, when atmospheric data are also included in the fit these 
ambiguities, although far from being resolved, become at least non-symmetric. In particular, 
the atmospheric bound on ^13 (blue lines) considerably depends on the mass hierarchy, and the 
global fit with atmospheric data included (red lines) exhibit a weak but visible dependence on 
the CP phase and on the octant. Of course, this is not a proof that atmospheric data will be 
relevant in the future, especially since the results of forthcoming accelerator experiments such as 
T2K will have a non-trivial dependence on the value of 5cp and on the neutrino mass hierarchy. 
However, it is a hint the very broad-range information provided by atmospheric neutrino data, 
which span about three orders of magnitude in length and more than five in energy, may still 
be complementary to accelerator experiments, which despite their high degree of precision are 
limited to a fixed baseline and cover only a very limited range in neutrino energy. In the rest of 
this talk we will present a systematic study of these potentialities. 

2 Sensitivity to oscillation parameters 

As already mentioned, the strength of atmospheric data is its very broad interval in neutrino 
energy {E^) and baseline (determined by the nadir angle &u)- In order to provide a global 
view of this whole range, in this section we will make extensive use of neutrino oscillograms of 




dependence on the neutrino mass hierarchy. 

the Earth, i.e. contours of equal probabihties in the neutrino energy-nadir angle plane P More 
exactly, we will show contours of equal oscillated-to-unoscillated ratio of atmospheric neutrino 
fluxes of flavor /3 arriving at the detector, $^/<I>Jj, obtained by folding the primary neutrino flux 
with the relevant conversion probability Pa^p, so that = Y^^^^ Pa^f^. As we will see, 
different regions of the {E,^, Q^) plane will show characteristic structures whose position and 
size is determined by various neutrino parameters. 

^13. Let us start by considering the sensitivity to ^13. As can be seen from Fig. |2|^a), for non- 
zero value of ^13 matter effects induce a resonance in the Vf. conversion probability at 
-Bi/ ~ 3 6 GeV. The precise position of this peak in the [Ev-, ^v) plane depends on the value 
of ^13, so that in principle atmospheric neutrinos could be used to measure this angle as long as 
it is larger than about 3°. However, in practice the sensitivity is limited by two factors: 

- statistics: at E^^ ^ 6 GeV the atmospheric flux is already considerably suppressed; 

- background: the i/^ — > Ve signal is diluted by the unavoidable background of I'e Vf. events. 
Therefore, although some sensitivity is to be expected in a megaton detector, it is likely that 
atmospheric neutrinos will not be competitive with dedicated long-baseline and reactor exper- 
iments for what concerns the determination of ^13. However, an explicit observation of this 
resonance will provide a very important confirmation of the MSW and parametric-resonance 
mechanisms. 

Hierarchy. As shown in Fig. [2]^b), the sensitivity to the hierarchy arises from the observation of 
the same high-energy resonance which is also involved in the sensitivity to ^13 . It is therefore only 
possible if ^13 is large enough. Note that in order to determine the hierarchy it is not sufficient 
to see the resonance (which would simply be a indication of non-zero ^13), it is also necessary to 
tell whether it occurred for neutrinos (normal hierarchy) or antineutrinos (inverted hierarchy). 
It is therefore particularly important to have a detector capable of charge discrimination. In 
the case of a non-magnetic detector such as a Water Cerenkov, if a resonance is observed it 
might still be possible to resolve the hierarchy by looking its size, since the number of neutrinos 
interacting in the detector is considerably larger than the number of antineutrinos and therefore 




a normal (inverted) hierarchy would result in a larger (smaller) signal. Note, however, that the 
amplitude of the peak is affected both by ^13 and by the value of ^23 , as we will see in the next 
paragraph, so that a poorly known value of these parameters will result in a considerable loss 
of sensitivity. 

Octant. The sensitivity to the octant is one of the topics where atmospheric neutrino are mostly 
useful. As can be seen Fig. [sj^a), we have two characteristic signatures: 

- at low energy {Ey < 1 GeV), we observe an excess (deficit) in the fe flux with respect to 
maximal mixing if $23 is smaller (larger) than 45° . This effect is due to subleading oscillations 
induced by Am^i, and is present also for ^13 = 0. For ^13 7^ the flux arriving at the 
detector is modulated with the very fast Am|^ oscillations, but the effect persists on average. 
Finally, this effect appears with the same sign for both neutrinos and antineutrinos, so that 
no charge discrimination is required for its identification. 

- at high energy {Ey > 3 GeV), we observe a decrease (increase) in the Ve flux with respect to 
maximal mixing if ^23 is smaller (larger) than 45°. This effect is again related to the matter 
resonance discussed for ^13 and the hierarchy, and indeed it appears only for ^13 7^ 0. As 
already seen, in a detector without charge discrimination the signal could be considerably 
suppressed. 

Note that the presence of a low-energy effect independent of ^13 guarantees a minimum sensitivity 
to the octant from atmospheric neutrinos, provided that the deviation of ^23 from maximal 
mixing is large enough. This is a unique feature which will prove very synergic with long- 
baseline data, as we will show in the next section. Note also that the slight preference for 
^23 < 45° visible in Fig. [T] arises precisely from this effect, and from the observation of a small 
excess in sub-GeV e-like events in Super-Kamiokande data. 

CP phase. Finally, let us spend a word on the sensitivity to the CP phase. A characteristic 
signal is visible in the intermediate-energy region, 1 GeV < E^ < 2> GeV, and arises from 
the interference of Am^i-induced and Am|]^-induced oscillations. Although in principle it is 
observable, as Fig. [T] demonstrates with present data, this effect is quite small and probably 
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Figure 4: Allowed regions in sin^ 26i3 and Sep for LBL data alone (contour lines) and LBL+ATM data com- 
bined (colored regions).'^ jjtr/wr Qtr/wr ].g£g]-g solutions with the true/wrong mass hierarchy and octant, 
respectively. The true parameter values are Sep — — O.SStt, sin'^ 612 — 0.3, sin^ 2^13 = 0.03, sin^ 623 = 0.6, 
Amii = 7.9 X 10"^ eV^ Amii = 2.4 x 10"=* eV^ 



hard to disentangle from other parameters. Moreover, its presence depends crucially on the size 
of 013, so as for the sensitivity to the hierarchy it is not "guaranteed". In the fits which we will 
discuss in the rest of this talk the impact of 6cp on the determination of the other parameters is 
properly taken into account, however we will not present any systematic study of the sensitivity 
of atmospheric data to 6cp itself. 



3 Synergies with long-baseline experiments 

So far we have discussed the potentialities of atmospheric neutrinos in general terms. Let us 
now consider concrete experimental setups and compare their performances. In particular, we 
will focus on three proposed experimentsP 

- a Beta Beam (/3B) from CERN to Frejus (130 Km). We assume 5 years of fg from ^^Ne and 
5 years of i>e from ^He at 7 = 100, with an average neutrino energy (E^) = 400 MeV. For the 
detector we assume the MEMPHYS Water-Cerenkov proposal, corresponding to 3 tanks of 
145 Kton each; 

- a Super Beam (SPL) from CERN to Frejus (130 Km). We assume 2 years of and 8 years of 
i'^ running, with an average energy {E,y) = 300 MeV. Again we use MEMPHYS as detector; 

- the T2K phase II (T2HK) experiment, corresponding to a 4MW super beam from Tokai to 
Kamioka (295 Km), with 2 years of and 8 years of z^^. The detector is the proposed 
Hyper-Kamiokande, rescaled to 440 Kton for a fair comparison with the /3B and the SPL. 

A characteristic feature in the analysis of future LBL experiments is the presence of parameter 
degeneracies. Due to the inherent three-flavor structure of the oscillation probabilities, for a given 
experiment in general several disconnected regions in the multi-dimensional space of oscillation 
parameters will be present. Traditionally these degeneracies are referred to as follows: 

- the intrinsic degeneracy: for a measurement based on the — > i^e oscillation probability for 
neutrinos and antineutrinos two disconnected solutions appear in the {Sep, 613) plane; 

- the hierarchy degeneracy: the two solutions corresponding to the two signs of Am|^ appear 
in general at different values of 6cp and ^13; 

- the octant degeneracy: since LBL experiments are sensitive mainly to sin^ 2^23 it is difficult 
to distinguish the two octants 623 < 45° and 823 > 45°. Again, the solutions corresponding 
to 023 and tt/2 — 623 appear in general at different values of 5cp and ^13. 



2o sensitivity to normal hierarchy from LBL + ATM data lbl+atm data to the octant of 
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Figure 5: Left: sensitivity to the mass hierarchy at 2a (Ax^ = 4) as a function of the true values of sin^ 20i3 and 
dcF- The solid curves are the sensitivities from the combination of long-baseline and atmospheric neutrino data, 
the dashed curves correspond to long-baseline data only. Right: Ax^ of the solution with the wrong octant of O23 
as a function of the true value of sin'^ 623- We have assumed a true value of ^13 — 0. 

This leads to an eight-fold ambiguity in ^13 and 5cp, and hence degeneracies provide a serious 
limitation for the determination of ^13, 5cp and the sign of Am|^. In Fig. |4] we illustrate 
this problem for the /3B, SPL and T2HK experiments. Assuming the true parameter values 
6cp = — 0.85-7r, sin^ 2^13 = 0.03 and sin^ ^23 = 0.6 we show the allowed regions in the plane 
of sin^ 2^13 and 6cp taking into account the solutions with the wrong hierarchy and the wrong 
octant. As visible in this figure, for the /3B the intrinsic degeneracy cannot be resolved, due 
to the poor spectral information and the lack of precise information on |Am|]^| and sin^ 2^23 
(usually provided by the t'^ disappearance), while for the super beam experiments SPL and 
T2HK there is only a four-fold degeneracy related to the sign of Am^^ and the octant of ^23- 
On the other hand, once atmospheric data are included in the fit all the degeneracies are nearly 
completely resolved, and the true solution is identified at 95% CL. This clearly show the presence 
of a synergy between atmospheric and long-baseline data: at least for this specific example, the 
combination of the two sets is much more powerful than the simple sum of each individual data 
sample. 

To further investigate this synergy, in the left panels of Fig. [5] we show how the combination of 
ATM-t-LBL data leads to a non-trivial sensitivity to the neutrino mass hierarchy. For LBL data 
alone (dashed curves) there is practically no sensitivity for the CERN-MEMPHYS experiments 
(because of the very small matter effects due to the relatively short baseline) , and the sensitivity 
of T2HK depends strongly on the true value of 6cp- However, by including the data from 
atmospheric neutrinos (solid curves) the mass hierarchy can be identified at 2a CL provided 
sin^ 2^13 > 0.02 0.03. As an example we have chosen in that figure a true value of ^23 = 
45°; in general the hierarchy sensitivity increases as ^23 increases .^1 Note that the sensitivity 
to the neutrino mass hierarchy shown in Fig. [5] is significantly improved with respect to our 
previous results.^l There are two main reasons for this better performance: first, we use now 
much more bins in charged lepton energy for fully contained single-ring events; second, we 
implemented also the information from multi-ring events. This latter point is important since 
the relative contribution of neutrinos and antineutrinos is different for single-ring and multi-ring 
events. Therefore, combining both data sets allows to obtain a discrimination between neutrino 
and antineutrino events on a statistical basis. This in turn contains crucial information on 
the hierarchy, since as discussed in Sec. [2] the mass hierarchy determines whether the matter 
enhancement occurs for neutrinos or for antineutrinos. 

In the right panel of Fig. [5] we show the potential of ATM-I-LBL data to exclude the octant 
degenerate solution. As seen in the previous section, this effect is based mainly on oscillations 
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Figure 6: Results for the reference fit for the angular averaged muon neutrino plus antineutrino flux extrapolated 
to the high-energy region compared to the corresponding data from AMANDA. Also shown are the results with 
different reference fluxes and the comparison with the predictions of different flux models. 



with Amg]^ and therefore we have very good sensitivity even for ^13 = 0; a finite value of 
^13 in general improves the sensitivity.^ From the figure one can read off that atmospheric 
data alone can resolve the correct octant at 3a if | sin^ ^23 — 0.5| > 0.085. If atmospheric 
data is combined with the LBL data from SPL or T2HK there is sensitivity to the octant for 
I sin'^ ^23 — 0.5 1 > 0.05. The improvement of the octant sensitivity with respect to previous 
analyse^30 follows from changes in the analysis of sub-GeV atmospheric events, where now 
three bins in lepton momentum are used instead of one. Note that since in this figure we have 
assumed a true value of ^13 = 0, combining the /3B with ATM does not improve the sensitivity 
with respect to atmospheric data alone. 



4 Direct determination of atmospheric fluxes 

So far we have discussed the potentialities of atmospheric neutrino for what concerns the de- 
termination of the neutrino parameters. The main message of the first part of this talk is that 
atmospheric data can provide very useful information on the neutrino mass matrix, despite the 
very large uncertainties in the neutrino fluxes. However, it is logically acceptable to invert the 
strategy, and to regard the poorly known atmospheric neutrino fluxes as a subject of investi- 
gation themselves. In this section we will therefore assume that the neutrino parameters have 
been accurately measured by other experiments, and we will show that it is possible to extract 
the the atmospheric neutrino fluxes directly from the data.l^ 

There are several motivations for such direct determination of the atmospheric neutrino 
fluxes. First of all it would provide a cross-check of the standard flux calculations as well 
as of the size of the associated uncertainties (which, being mostly theoretical, are difficult to 
quantify). Second, a precise knowledge of the atmospheric neutrino flux is of importance for 
high energy neutrino telescopes, both because they are the main background and they are used 
for detector calibration. Finally, such program may quantitatively expand the physics potential 
of future atmospheric neutrino experiments. Technically, however, this program is challenged 
by the absence of a generic parametrization of the energy and angular functional dependence of 
the fluxes which is valid in all the range of energies where there is available data. We bypass this 
problem by using artificial neural networks as unbiased interpolants for the unknown neutrino 
fluxes. However, the precision of the available experimental data is not yet enough to allow for 



a separate determination of the energy, zenith angle and flavor dependence of the atmospheric 
flux. Consequently in our work we have assumed the zenith and flavor dependence of the flux 
to be known with some precision and extract from the data only its energy dependence. Thus 
the neural network flux parametrization will be: 

^1]^{E,, ca, h) = F'^'iE,) <^':%{E,, Ca, h) (4) 

where F^^{Eiy) is the neural network output when the input is the neutrino energy E^y. 

In Fig. |6]we show the results of our flt to the atmospheric neutrino flux as compared with the 
computations of the HondsPand BartoP groups. The results of the neural network flt are shown 
in the form of a Icr band, and plotted as a function of the neutrino energy. For comparison we 
also show the data from AMANDA.!^ We see from this flgure that the flux obtained from the 
flt is in reasonable agreement with the theoretical calculations, although the flt seems to prefer 
a slightly higher flux at higher energies. This indicates that until about E^ ^ \ TeV we have a 
good understanding of the normalization of the fluxes, and that the present accuracy from Super- 
Kamiokande neutrino data is comparable with the theoretical uncertainties from the numerical 
calculations. Note also that the results of our alternative flts depends only mildly on the choice 
of Honda or Bartol as the reference flux. This suggests that the present uncertainties on the 
angular dependence have been properly estimated, so that the assumed angular dependence 
has very little effect on the determination of the energy dependence of the fluxes. Thus the 
atmospheric neutrino flux determined with our method could be used as an alternative to the 
existing flux calculations. 



5 Conclusions 



In this talk we have discussed the potentialities of atmospheric neutrino data in the context of 
future neutrino experiments. We have shown that despite the large uncertainties in the neutrino 
fluxes atmospheric data will still provide useful information on the neutrino parameters, due 
to their very broad range in neutrino energy and nadir angle. In particular, we have proved 
that the sensitivity obtained by a combination of atmospheric and long-baseline data is much 
stronger than the one achievable by each data set separately. Finally, we have shown that present 
atmospheric data can be used to obtain a direct determination of the atmospheric neutrino fluxes. 
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